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Abstract 



Physics of radioactive nuclear beams is one of the main frontiers of nuclear 
science today. Experimentally, thanks to technological developments, we are 
Q\ \ on the verge of invading the territory of extreme N/Z ratios in an unprece- 

dented way. Theoretically, nuclear exotica represent a formidable challenge 
■ for the nuclear many-body theories and their power to predict nuclear prop- 

erties in nuclear terra incognita. It is important to remember that the lesson 
\ learned by going to the limits of the nuclear binding is also important for 

"normal" nuclei from the neighborhood of the beta stability valley. And, of 
course, radioactive nuclei are crucial astrophysically; they pave the highway 
"q ■ along which the nuclear material is transported up in the proton and neutron 

numbers during the complicated synthesis process in stars. 



I. INTRODUCTION 

There are only 263 stable nuclei; they are surrounded by radioactive ones. Some of the 
unstable nuclei are long-lived and can be found on Earth, some are man-made (actually, as 
many as ~2,200 nuclei have been produced in nuclear laboratories), and several thousand 
nuclei are the yet-unexplored exotic species. The decay characteristics of most radioactive 
nuclei are determined by the beta decay, i.e., by weak interactions. For heavier nuclei, where 
the electromagnetic interaction plays a more important role, other decay channels, such as 
emission of alpha particles or spontaneous fission, dominate. Moving away from stable nuclei 
by adding either protons or neutrons, one finally reaches the particle drip lines. The nuclei 
beyond the drip lines are unbound to nucleon emission; that is, for those systems the strong 
interaction is unable to bind A nucleons as one nucleus. 



* Will be published in Theme Issue on Science with Beams of Radioactive Nuclei, Philosophical 
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So, the territory of exotic nuclei is enormous. The uncharted regions of the (N, Z) plane 
contain information that can answer many questions of fundamental importance for science: 
How many protons and neutrons can be clustered together by the strong interaction to form a 
bound nucleus? What are the proton and neutron magic numbers of the exotic nuclei? What 
are the properties of very short-lived exotic nuclei with extreme neutron-to-proton ratio 
N/Z7 What is the effective nucleon-nucleon interaction in the nucleus having a very large 
neutron excess? There are also related questions in the field of nuclear astrophysics. Since 
radioactive nuclei are produced in many astrophysical sites, knowledge of their properties is 
crucial to the understanding of the underlying processes. 

Nuclear life far from stability is different from that around the stability line; the promised 
access to completely new combinations of proton and neutron numbers offers prospects 
for new structural phenomena. The unique structural factor is the weak binding; hence 
closeness to the particle continuum. The main objective of this paper is to discuss some of 
the theoretical challenges and opportunities of research with exotic nuclear beams. 

The paper is organized as follows. Section [II] contains a brief review of theoretical de- 
velopments related to the physics of exotic nuclei. Sections [ITT] and [TV] contain some of the 
physics issues of the neutron and proton drip lines, respectively, including a discussion on 
the theoretical uncertainties in determining the particle drip lines. Finally, conclusions are 
contained in Sec. 0. 



II. NEW THEORETICAL ASPECTS OF PHYSICS WITH EXOTIC BEAMS 



From a theoretical point of view, spectroscopy of exotic nuclei offers a unique test of 
those components of effective interactions that depend on the isospin degrees of freedom. 
In principle, the effective nucleon-nucleon interaction in heavy nuclei should be obtained by 
means of the Bruckner renormalization which corrects the free interaction for the effects due 
to the nuclear medium. In practice, however, the effective interaction is approximated by 
means of some phenomenological density-dependent force with parameters that are usually 
fitted to stable nuclei and to selected properties of the infinite nuclear matter. Hence, it 
is by no means obvious that the isotopic trends far from stability predicted by commonly 
used effective interactions are correct. In the models aiming at such an extrapolation, the 
important questions asked are: What is the density dependence of the two-body central 
force [|l]-||? What is the density and radial dependence of the one-body spin-orbit force 
|Q-0]? Does the spin-orbit splitting strongly vary with N/Z H? What is the form of the 
pairing interaction in weakly bound nuclei |9Hl2|? What is the importance of the effective 
mass (i.e., the non- locality of the force) for isotopic trends [0? What is the role of the 
medium effects (renormalization) and of the core polarization in the nuclear exterior (halo 



or skin region) where the nucleonic density is small ||14|| ? Similar questions are asked in 



connection with properties of nuclear matter [|6||l5|-|T7|j , neutron droplets pfl, and the physics 



of the neutron-star crust 18,19 



The radioactive nuclear beams experimentation is expected to expand the range of nuclei 
known. That is, by going to nuclei with extreme N/Z ratios, one can magnify the isospin- 
dependent terms of the effective interaction (which are small in "normal" nuclei). The hope 
is that after probing these terms at the limits of extreme isospin, we can later go back to 
the valley of stability and improve the description of "normal" nuclei. 
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But this task is not going to be easy. In many respects, weakly bound nuclei are much 
more difficult to treat theoretically than well-bound systems. Hence, before tackling the 
problem of force parametrization at the extremes, one should be sure that the applied 
theoretical tools of the nuclear many-body problem are appropriate. 

As mentioned above, the main theoretical challenge is the correct treatment of the par- 
ticle continuum. For weakly bound nuclei, the Fermi energy lies very close to zero, and 
the decay channels must be taken into account explicitly. As a result, many cherished 
approaches of nuclear theory such as the conventional shell model, the pairing theory, or 
the macroscopic-microscopic approach must be modified. But there is also a splendid op- 
portunity: the explicit coupling between bound states and continuum, and the presence 
of low-lying scattering states invite strong interplay and cross-fertilization between nuclear 
structure and reaction theory. Many methods developed by reaction theory can now be 
applied to structure aspects of loosely bound systems. 

How to extend traditional tools of nuclear theory to account for the scattering of nucleons 
from bound single-particle orbitals to unbound states? The closeness of particle continuum 
reverberates in two aspects of the theoretical description. Firstly, the particles forming a 
bound nuclear state can effortlessly virtually scatter back and forth into the particle contin- 
uum phase space. This process must conserve the compactness of the nuclear wave function 
which remains bound even with such a virtual scattering taken into account. A theoretical 
description of this kind of effects still remains a virgin territory, although some progress 
has been made in the analysis of the virtual pair scattering pO|JlQ[| . Secondly, nucleons can 



very easily leave the nucleus altogether and enter the particle continuum through the real 
scattering. For this, it is enough to slightly shake the nucleus by providing it with a little bit 
of energy. This is an old problem which, in the context of excited states near or above the 
particle threshold, has been addressed by the continuum shell model (CSM) [f2~T|-|2"E|] . In the 
CSM, the continuum states (decay channels) and bound states are treated on equal footing. 
Consequently, correlations due to the coupling to resonances, the spatial extension effects 
in weakly bound states, the structure of resonances, and the structure of particle transfer 
form factors are properly described by the CSM. So far, most applications of the CSM have 
been concerned with the situation when there is only one particle occupying the shell-model 



continuum. This is because the continuum-continuum coupling is difficult to treat p9 



Often, particle continuum is approximated by the quasibound states, i.e., the states re- 
sulting from the diagonalization of a finite potential in a large basis |30| , |3T|| or by enclosing 
the finite nuclear potential within an infinite well with walls positioned at a large distance 
from the nuclear surface [p2| , p3[ . More sophisticated methods of discretizing continuum in- 
clude the Sturmian function expansions and resonant state expansions. Sturmian functions, 
also known as Weinberg states, form a discrete set of states which behave asymptotically 
as outgoing waves. They have been used as a basis in the solution of scattering equations, 



including various applications of the CSM [22,33-R^]. The Gamow states are eigenstates of 



the time-independent Schrodinger equation with complex eigenvalues |38|-f|(J. They have 
been applied to many problems involving an unbound spectrum [¥I-i3"|. 

In the description of weakly bound systems, pairing interaction plays a unique role (see 
Section [111B| below). In the BCS or BCS-like methods based on bound and quasibound 



states (for examples of such calculations see, e.g., Refs. [[H],[|4[]), the virtual scattering of 



nucleonic pairs from bound states to the positive-energy states leads to the presence of 
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a "particle gas" surrounding the nucleus fl20f . To show it, we decompose the BCS wave 
function into contributions from bound states (e< < 0) and quasibound states (e 3 - > 0): 



^bcs= n 

i,e;<0 



Us + ViaSat 



n 

3,tj>0 



Ui + V ja }a\) |0> 



(1) 



While the bound-state component in Eq. (|T]) represents the localized wave function, i.e., it 
decays asymptotically, the second part represents the contribution from quasibound states 
and leads to non-localized densities with incorrect asymptotic behavior. Indeed, although 
the nuclear densities eventually vanish at large distances by construction (finite size of the 
basis, finite size of the box in which calculations are performed), the wave functions of 
positive-energy states do not decay outside the nuclear volume. As discussed below in 
Sec. 1 1 j~Bl this problem is overcome in the Hartree-Fock-Bogolyubov (HFB) method with a 
realistic pairing interaction in which the coupling of bound states to the particle continuum 
is correctly taken into account |20|JT0| . 
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FIG. 1. Top: Two- neutron separation energies for the even-even nickel isotopes predicted in 
the HFB and HF+BCS calculations with the Skyrme interaction SkP. Bottom: The predicted rms 
neutron radii. The arrows indicate the neutron subshell closures at N=28 and 50. (From Ref. 



To illustrate this point, critical in the context of calculations for drip- line nuclei, Fig. [I 
displays theoretical two-neutron separation energies, S^n, and the rms neutron radii for the 
even-even nickel isotopes obtained in the HFB and HF+BCS calculations. In the HF+BCS 
variant, the self-consistent pairing gaps obtained from the HFB calculations were used within 
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the fixed-gap approximation. As seen in Fig. |T|, the values of S2 n obtained in the HFB and 
HF+BCS calculations agree very well; some deviations are seen only for the neutron drip- line 
systems with A>60, where the HF+BCS values are slightly lower. However, this excellent 
agreement does not extend to neutron radii. In the HFB calculations the neutron radii 
behave very smoothly as a function of N. On the other hand, in the HF+BCS model there 
is a dramatic increase in the neutron radii between magic numbers and for the neutron-rich 
nuclei resulting from the unphysical occupation of positive-energy quasibound states. The 
difference between values of radii obtained within the HFB and HF+BCS calculations can 
be as large as 0.8 fm and increases dramatically for weakly bound systems. Other examples 



illustrating the unphysical effect of the particle gas can be found in Refs. |10| , |^6| . 

Consequently, for large exotic nuclei, the self-consistent HFB treatment is not a matter of 
choice, it is a must. The calculations are not easy, especially if the self-consistent symmetries 
(e.g., spherical symmetry) are broken. Possible strategies for solving the HFB equations in- 
clude the two-step diagonalization ||47|| , the gradient method in the canonical representation 



48] , or the state-dependent Hamiltonian method [ 49fl . Some examples of HFB calculations 



are discussed in the following sections. Other calculations with self-consistent inclusion 
of pairing and continuum states include: quasiclassical Lagrangian method calculations in 
the coordinate representation [^DJ-Q and application of the relativistic Hartree-Bogolyubov 
theory in coordinate space to light nuclei [55|.56]. 

In order to describe excited states, one has to go beyond the mean-field approximation. 
One method, often used in the description of low-lying collective states in drip-line nuclei, 
is the continuum random phase approximation (CRPA) based on the single-particle Green 
function approach in the coordinate representation |57]J58| . In particular, the CRPA has been 



employed extensively to the very low-energy multipole strength in drip-line nuclei [59-62 



However, for meaningful predictions of excited states in weakly bound nuclei where pairing 
is expected to be very important, it is necessary to use the quasiparticle RPA scheme based 
on the coordinate space HFB formalism. Calculations along these lines can be found in 



recent Refs. |63],|64]. 



The following sections contain a brief description of selected highlights of the physics of 
radioactive nuclear beams. It is to be noted that there are many other important topics 
that have been left out in our discussion (e.g., physics of very light drip-line nuclei, physics 
of beta decay, reaction aspects involving radioactive ions) . A key point is that the variety 
of exciting new phenomena is one of the driving forces behind research with exotic beams. 



III. PHYSICS OF NEUTRON-RICH NUCLEI: UPPER LIMITS OF THE N/Z 

RATIO 

Since neutrons do not carry an electric charge and do not repel each other, many neutrons 
can be added to nuclei starting from the valley of stability. As a result, the "lever arm" 
separating the neutron drip line from the valley of stability is large and difficult to probe 
experimentally; except for the lightest nuclei, the bounds of the neutron stability are not 
known. But it is just for the nuclei with the extreme neutron excess that theory predicts 
many new and initially unexpected phenomena. In addition to nuclear structure interest, 
the neutron-rich environment is important for astrophysics and cosmology. 
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A. Density distributions 



Neutron-rich nuclei are characterized by spatially extended density distributions which 
give rise to large nuclear radial moments. Extreme cases are halo nuclei - loosely bound 
few-body systems with about thrice more neutrons than protons. The halo region is a zone 
of weak binding in which quantum effects play a critical role in distributing nuclear density 
in regions not classically allowed. 

In the heavier, neutron-rich nuclei, where the concept of mean field is better applicable, 
the separation into a "core" and "valence nucleons" seems less justified. However, also 
in these nuclei the weak neutron binding implies the existence of the neutron skin (i.e., a 
dramatic excess of neutrons at large distances). 
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FIG. 2. Self-consistent spherical neutron densities Pn(t) calculated with the SkP, SIII* 5 , and 
D1S interactions for selected tin isotopes across the /3-stability valley. Since proton and neutron 
densities in the nucleus 100 Sn are very similar p5f , the deviation from the neutron density at iV=50 
roughly represents the skin effect. (From Ref. 

Figure 



displays the neutron HFB densities for several tin isotopes across the stability 

Slir 5 0, and D IS [| 



valley, calculated with the effective interactions SkP |20 



]. The 

densities obtained with these forces are qualitatively very similar. One can see that adding 
neutrons results in a simultaneous increase of the central neutron density, and of the density 
in the surface region. The relative magnitude of the two effects is governed by a balance 
between the volume and the surface asymmetry energies of effective interactions. Since all 
three forces considered have been fitted in a similar way to bulk nuclear properties (including 
the isospin dependence), the resulting balance between the volume and the surface isospin 
effects is similar. Of course, this does not exclude some differences which are seen when a 
more detailed comparison is carried out. As will be seen in Sec [T I [TIJ below, pair densities 
depend much stronger on the effective forces. 
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FIG. 3. Comparison of experimental (symbols) and calculated (HFB+SkP, solid line; 
HFB+SIir^, dashed line) proton and neutron rms radii for the semi-magic nuclei with N (or 
Z)=20, 28, 50, and 82. The experimental proton radii were extracted from experimental charge 



radii [S7|. The experimental neutron radii were obtained from the analysis of the high-energy 
proton- nucleus scattering ||68[|. (From Ref. p5|.) 



A very interesting aspect of nuclei far from stability is an increase in their radial di- 
mension with decreasing particle separation energy |69| . This effect is especially strong in 
weakly bound nuclei close to the neutron drip line. Figure || presents the comparison of 
experimental (symbols) and calculated proton and neutron rms radii for the semi-magic 
isotopes and isotones with Z (or N)=20, 28, 50, and 82. For neutron radii in general, and 
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for proton radii in light and medium-mass nuclei, the HFB+SkP and HFB+SIII^ models 
yield very similar predictions. In heavy nuclei (Z or iV=50, 82) the proton radii in SkP are 
lower than in SIII 5p , in closer agreement with experiment. In Fig. |3], the neutron skin is 
manifested by a rapid increase in neutron rms radii when approaching the two-neutron drip 
line. This effect is very localized in neutron number; it appears only in a few nuclei in the 
immediate neighborhood of the neutron drip line. 



B. Pairing correlations 



Pairing correlations play a very special role in drip-line nuclei [20|,10|. This is seen from 
the approximate HFB relation between the Fermi level A, pairing gap A, and the particle 
separation energy S ~ —A — A. At the drip line 5* is very small and A + A ~ 0. Consequently, 
the single-particle field characterized by A and the pairing field A are equally important. 
In other words, contrary to the situation encountered close to the line of beta stability, the 
pairing component of the Hamiltonian can no longer be treated as a residual interaction, 
i.e., a small perturbation important only in the neighborhood of the Fermi surface. 

Surprisingly, very little is known about the pairing channel of the effective nucleon- 
nucleon interaction. In most calculations, the pairing Hamiltonian has been approximated 
by the state-independent seniority pairing force, or schematic multipole pairing interaction 
|70| . Such oversimplified forces, usually treated by means of the BCS approximation, perform 



remarkably well when applied to nuclei in the neighborhood of the stability valley (where, as 
pointed out above, pairing can be considered as a small correction). As a result, considerable 
effort was devoted in the past to optimizing the HF part of the interaction, while leaving 
the pairing component aside. 

A detailed discussion of the present status of effective interactions in the particle-particle 



channel can be found in Ref. [10|. The main questions pertaining to this problem are: What 



is the microscopic origin of the pairing interaction ||7lH75||? What is the role of finite range 



and the importance of density dependence |76|-pT|? How can properties of the pairing force 



be tested experimentally? These questions are of considerable importance not only for 
nuclear physics but also for nuclear astrophysics and cosmology |18],^2]|33] . 



Because of strong surface effects, the properties of weakly bound nuclei are perfect lab- 
oratories in which to study the density dependence of pairing interactions. As an example 
of what can be expected far from stability, Fig. |] displays the neutron HFB local pairing 
densities p(r) calculated for several tin isotopes across the stability valley, and for three dif- 
ferent effective interactions: Skyrme interactions SkP and SkP* 5 , and the finite-range Gogny 
interaction D1S. [The density p(r) is proportional to the probability of finding the correlated 
pair of nucleons at point r, see Ref. [|10| for definitions and discussion.] The pairing densi- 



ties shown in Fig. ^ nicely reflect different characters of the interactions used. Namely, the 
contact force SkP 5 leads to pairing densities that are, in general, largest at the origin and 
decrease towards the surface; this is characteristic of the volume-type pairing correlations. 
A different pattern appears for the SkP results, where the density dependence renders the 
pairing interaction strongly peaked at the surface. In this case, the pairing densities tend to 
increase when going from the center of the nucleus towards its surface. A more pronounced 
dependence on the neutron excess is seen here in the surface region. Near the drip line, 
the pairing density develops a long tail extending towards large distances. The results ob- 



8 



tained for the finite-range Gogny interaction exhibit features intermediate between surface 
and volume-type pairing correlations. 
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FIG. 4. Self-consistent spherical neutron pairing densities /5j\r(r) calculated with the SkP, 
SkP 5 , and D1S interactions for selected tin isotopes across the /3-stability valley. (From Ref. ]i0[.) 

An experimental observable that may probe the character of the pairing field is the 
pair transfer form factor, directly related to the pairing density p. The difference in the 
asymptotic behavior of single-particle density p and pair density p in a weakly bound system 
can be probed by comparing the energy dependence of one-particle and pair-transfer cross 
sections. Such measurements, when performed for both stable and neutron-rich nuclei, can 
shed some light on the asymptotic properties of the pair densities; hence on the character 
of the pairing field. 

Figure [| displays the pair transfer form factors r 2 p(r) calculated in 120 Sn, 150 Sn, and 
172 Sn with the SkP interaction. These pair transfer form factors clearly show that this 
process has a predominantly surface character. In particular, there is a significant increase 
in the pair transfer form factors in the outer regions of drip-line nuclei. In 120 Sn, the form 
factors vanish around 9 fm, while in 150 Sn and 172 Sn they extend to much larger distances. 



C. Shell structure far from stability and position of the neutron drip line 

The structure of nuclei is expected to change significantly as the limit of nuclear stability 
is approached in neutron excess. Due to the systematic variation in the spatial distribution 
of nucleonic densities and the increased importance of the pairing field, the average nucleonic 
potential is modified when approaching the neutron drip line. The main effect is the increase 
of the potential diffuseness; the single-particle neutron potential in drip-line nuclei becomes 
very shallow, and the resulting shell-model spectrum resembles that of a harmonic oscillator 
with a spin-orbit term and with a weakened £ 2 term [65]. This results in a new shell structure 



characterized by a more uniform distribution of normal-parity orbits and the unique-parity 
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intruder orbit which reverts towards its parent shell. 
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FIG. 5. Pair transfer form factor, r 2 p(r), calculated directly from the HFB pairing density 
p(r). (From Ref. p 



The effect of the weakening of shell effects in drip-line nuclei, first mentioned in Ref. 
was further investigated in Refs. |85|j86| , |65| . 



8|. Quenching of shell effects manifests 
itself in the behavior of two-neutron separation energies S^. This is illustrated in Fig. |6| 
which displays the two-neutron separation energies for the iV=80, 82, 84, and 86 spherical 
even-even isotones calculated in the HFB model with the SkP [20] and SLy4 [^{J effective 
interactions. The large N=82 magic gap, clearly seen in the nuclei close to the stability 
valley and to the proton drip line, gradually closes down when approaching the neutron 
drip line. This result is independent of the size of the N=82 shell gap in stable nuclei, 
which is slightly underestimated and overestimated by SkP and SLy4 forces, respectively, as 
compared to experimental data. It can be attributed to two effects: (i) a gradual increase 
of the neutron surface diffuseness across the stability valley related to an increase of the 
neutron excess, and (ii) the influence of the continuum, which results in closing the shell 
gap near the neutron drip line down to zero. 
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FIG. 6. Two-neutron separation energies for the N=80, 82, 84, and 86 spherical even-even 
isotones calculated in the HFB+SkP and HFB+SLy4 <5 models as functions of the proton number. 
The arrows indicate the proximity of neutron and proton drip lines for small and large proton 
numbers, respectively. 

Predicted behaviour of the two-neutron separation energies depends very much on the 
effective interaction used. This is illustrated in Fig. [7] which shows the Sin values calculated 
in tin isotopes for the Gogny interaction D1S and for four variants of the Skyrme interaction. 
The Gogny force and the SkP and SLy4 <5 Skyrme forces predict a gradual decrease of the 
two- neutron separation energies while the older Skyrme forces, SIII 5 and SkM* 5 , give almost 
constant values followed by a sudden drop at iV = 126. This shows that the iV=126 shell 
quenching is not a generic effect. 

As seen in Fig. [7], the position of the neutron drip line for the Sn isotopes also de- 
pends on the effective interaction used; it varies between iV=118 (D1S) and iV=126 (SIII* 5 ). 
Hence, even if the theoretical method used to calculate nuclear masses is reliable near the 
drip line (this is not the case for the commonly used macroscopic-microscopic models, see 
Ref. [[HJ]), the uncertainty due to the largely unknown isospin dependence of the force gives 
an appreciable theoretical "error bar" . 
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FIG. 7. Two-neutron separation energies S^n for the Sn isotopes, calculated in the HFB 
approach with the Skyrme interactions SIII* 5 , SkM 5 , SLy4 5 , and SkP and with the Gogny-DlS 
interaction. 



Unfortunately, the results presented in Fig. [7| do not tell us much about which of the forces 
discussed should be the "preferred one" since we are dealing with dramatic extrapolations 
far beyond the region known experimentally. The comparison with the data can be carried 
out for nuclei closer to stability, see Ref. |JD| for a recent quantitative analysis. A comparison 
with the experimental two-neutron separation energies is displayed in Fig. |[ As seen, the 
agreement with data is unsatisfactory for SIII 5 . In particular, the shell-gap sizes at N=50, 
82, and 126 are strongly overestimated, and the values and the slopes of S2 n are, in most 
cases, incorrect. For the SkM 5 force one obtains similar deficiencies ||13|| . The results shown 



in Figs. [7] and |8| illustrate a very strong dependence of the two-neutron separation energies 
on the force parameters. Although older forces, such as SIII and SkM*, can perform well in 
certain regions of Z and N, they do not give a satisfactory global reproduction of the data. 
On the other hand, a fairly good global agreement obtained with the SkP and SLy4' 5 suggests 
that the improvement is possible, while a still better parametrization would be welcome. Of 
course, forces which fail in reproducing the behavior with (N — Z) in known nuclei have little 
chance to perform better when going far from stability. For example, the predicted values 
of S2n obtained for the SIII* 5 and SkM* 5 interactions (Fig. [F]) followed by a strong shell effect 
at iV=126 do not seem very reliable. A detailed analysis of the force-dependence of results 
may give us valuable information on the relative importance of various force parameters. 

The gradual change in shell structure is expected to give rise to new sorts of collective 
phenomena pl|j9T[] . It is also to be noted that the experimentally observed collapse of magic 
gaps seen in some neutron-rich light nuclei is conventionally explained in terms of the shape 
transition to the deformed intruder configuration. Here, spectacular examples are 32 Mg 2 o 
92| , |93H and 44 S 2 s P^| , p5| . In both cases, HF calculations predict the shape transition |96|j97 



due to the crossing of spherical configuration by intruder states. Such a lowering of the in- 
truder configuration depends on the detailed balance between three components in the total 
energy: (i) the position of the intruder state at the spherical shape, (ii) the deformation en- 
ergy gain associated with the deformation-driving orbital, and (hi) the symmetry-restoring 
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force exerted by particles in the magic-shell configuration. The first and the third of these 
elements directly depend on the spherical shell structure and thus can be affected by the 
hypothetical shell quenching mechanism. However, the appearance of a coexisting configu- 
ration cannot, of course, be discussed solely in terms of the spherical shell structure, and it 
still remains an open and exciting problem. 
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FIG. 8. Two-neutron separation energies, S^n, for proton- magic isotopes. The HFB results 
with the SIII 5 , SLy^, and SkP parametrizations (solid lines) are compared with the experimental 
data (full circles) and systematic trends (open circles) |98|j . 

In the following section we discuss another important aspect of the shell quenching, 
namely, its consequences for the r-process and the stellar nucleosynthesis. 
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D. Structure of neutron-rich nuclei and the r-process 



The very neutron-rich drip-line nuclei cannot be reached experimentally under present 
laboratory conditions. On the other hand, these systems are the building blocks of the 
astrophysical r-process; their separation energies, decay rates, and neutron capture cross 
sections are the basic quantities determining the results of nuclear reaction network cal- 
culations. Consequently, one hopes to learn about properties of very neutron-rich systems 
by studying the r-process component of the solar-system abundances of heavy elements 
|99| , |100| , |88| , |101| , |102|| . The recent r-process network calculations [p^ , ^ , |lOT[] , based on several 
mass formulae, indicate that a quenching of the shell effect at iV=82 is required in order to 
fill the v4=120 and 140 r-abundance troughs, in accordance with the results of the HFB+SkP 
model shown in Fig. || 

In addition to nuclear structure, there are other factors which can influence the r-process 
abundances, for instance the astrophysical conditions of temperature, neutron density, and 
the process time scale ||103|| . The possibility that abundances of r-process elements may be 



altered by the intense neutrino flux has been discussed in Refs. [|104j , [105|| . According to their 
calculations, neutrino reactions can be important in breaking through the waiting-point 
nuclei at N=50 and 82, and the r-process abundances in the A=125 and 185 regions can be 
affected by neutrino post-processing effects. 

In addition to nuclear masses, another important piece of nuclear structure that deter- 
mines the path of the r-process is the GT strength. It is usually calculated in the quasiparticle 
RPA (QRPA) theory p9| , |106| - [108|| . The important development, yet to be done, would con- 
sist of performing systematic microscopic QRPA calculations based on the HFB densities. 
This would guarantee the proper treatment of the particle continuum in the weakly bound 
nuclei on the r-process path. 



E. Deformation of drip-line nuclei 



as well 



Neutron halos and heavy, weakly bound neutron-rich nuclei offer an opportunity to study 
the wealth of phenomena associated with the closeness of the particle threshold: particle 
emission (ionization to the continuum) and characteristic behavior of cross sections 
existence of soft collective modes and low-lying transition strength [|1 10| — 11 1 3| . |6Tf 
as various other nuclear properties in the sub-threshold regime. We have learned that 
weakly bound nuclei are different; they have giant sizes, they are diffused, they are strongly 
superfluid, their shell structure is probably different. But can they be deformed? 

The importance of non-spherical intrinsic shapes in halo nuclei has been discussed, es- 
pecially in the context of a one-neutron halo in n Be. The ground state of n Be is a l/2 + 
state. The low neutron separation energy, S' n =504 keV, allows for only one bound excited 
level (l/2~ at 320 keV). The halo character of n Be has been confirmed by studies of reac- 
tion cross sections ||1 14|| , and the importance of deformation can be inferred from the large 
quadrupole moment of its core 10 Be, |Q|=229mb [115]. The role of deformation in lowering 
the excitation energy of the l/2 + intruder level in n Be has been recognized ||116|J117|| , but 
the joint effect of loose binding and deformation has not been considered. (See, however, 
recent references [|118| - |122|j .) 
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In a recent study [ |123|| , the notion of shape deformations in halo nuclei has been addressed 
by considering the single-particle motion in the axial spheroidal square well. The properties 
of the deformed single-particle states, especially in the subthreshold region, were analyzed by 
making the multipole decomposition in the spherical partial waves with well-defined orbital 
angular momentum. It has been concluded that in the limit of very weak binding, the 
geometric interpretation of shape deformation is lost. That is, the deformation of the halo is 
solely determined by the spatial structure of the valence state wave function, independently 
of the shape of the core. The deformed core merely establishes the quantization axis of the 
system - important for determining the angular momentum projection on the symmetry 
axis, A. 

Figure || shows the contour map of Pg (probability to find the partial wave £ in a given 
Nilsson state [n exc A7r]) for the A=0 orbitals as functions of binding energy and deforma- 
tion. The structure of the [10+] Nilsson level, originating from the spherical Is state, is 
completely dominated by the £=0 component, even at very large deformations. A rather 
interesting pattern is seen in the diagram for the [20+] orbital originating from the spher- 
ical Id state. The £=2 component dominates at low and medium deformations, and the 
corresponding probability Pe=2 slowly decreases with 5 at large deformations approaching 
the (constant) asymptotic limit. However, a similar effect, namely the decrease of the £=2 
component, is seen when approaching the zero binding energy threshold. In the language of 
the perturbation theory ||21|| , this rapid transition comes from the coupling to the low-energy 
£=0 continuum. As a consequence of the dominating role of the s-wave (and p-wave, for 
negative-parity states), in the limit of weak binding, the total quadrupole deformation of 
the (core+valence) system depends solely on the geometry of the valence orbital. Namely, 
it is consistent with a superdeformed shape (vr=-, A=0 halo), a spherical shape (vr=+, A=0 
halo), or an oblate shape (tt=-, A=l halo), regardless of the deformation of the core. In 
the language of the self-consistent mean-field theory, this result reflects the extreme softness 
of the system to the quadrupole distortion. Shape deformation is an extremely powerful 
concept provided that the nuclear surface can be properly defined. However, for very dif- 
fused and spatially extended systems, the geometric interpretation of multipole moments 
and deformations is lost. 

The presence of the spatially extended neutron halo gives rise to the low-energy isovector 
modes. The deformation decoupling of the halo implies that the nuclei close to the neu- 
tron drip line are excellent candidates for isovector quadrupole deformations, with different 
quadrupole deformations for protons and neutrons. An example of such a situation has been 
predicted in the self-consistent calculations for the neutron-rich sulfur isotopes performed 
using the Skyrme HF and relativistic mean-field methods [ P7| , |124j| . When approaching the 
neutron drip line, the calculated values of quadrupole deformation for neutrons are system- 
atically smaller than those of the proton distribution. Another promising candidate for such 
effects is the 120 Sr nucleus which is presently being studied by the GCM configuration mixing 



of the Skyrme HF states [125 



The discussion of deformation in halo nuclei applies to systems with very small binding 
energy and with negligible pairing. How will this scenario be modified in the presence of 
pairing, and for greater separation energies? We do not know the answer to this question 
at present; the systematic investigation of the interplay between deformation and pairing in 
weakly bound neutron-rich nuclei is one of the most exciting avenues of RNB physics. 
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FIG. 9. Contour maps of probabilities Po an< i Pi f° r the [10+] (top) and [20+] (bottom) 
Nilsson levels as functions of deformation and binding energy. (From Ref. [p.23|.) 



IV. PHYSICS OF PROTON-RICH NUCLEI: AT AND BELOW THE N/Z=l 

LIMIT 

On the proton-rich side of the valley of stability, physics is different than in nuclei with 
a large neutron excess. Because of the Coulomb barrier which tends to localize the proton 
density in the nuclear interior, nuclei beyond the proton drip line are quasibound with respect 
to proton emission. However, in spite of the stabilizing effect of the Coulomb barrier, the 
effects associated with the weak binding are also present in proton drip-line nuclei. They 
are not as dramatic as on the other side of the stability valley, but nevertheless important. 
For instance, the Thomas- Ehrman shift [ 126 , |127|| , which is due to changes in the Coulomb 



energy of the weakly bound proton, can lead to a decrease in the energy differences between 
analog states by a few hundred keV. This effect is most significant for loosely bound states 



and for orbitals having low angular momentum | 128 |. (For the self-consistent calculations 
of the Thomas- Ehrman shift in the doubly-magic diproton emitter 48 Ni, see Ref. [|46|].) 
Consequently, indiscriminate applications of the nuclear shell-model to nuclei close to and 
beyond the proton drip line, ignoring the systematic changes in single-particle energies and 
wave functions due to weak binding, should probably be taken with a grain of salt. 

The doubly magic N=Z=50 nucleus 100 Sn is a paradigm of RNB physics at the proton- 
rich side. Although it was found experimentally three years ago [ |129|J130|| , it took more than 
two years to roughly determine its mass [|131|] , and it will probably take quite a few years to 
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find its first excited state. Actually, the question what is this state constitutes an unresolved 
problem which is a challenge for theoretical predictions. 



A. Neutron-proton correlations 



A unique aspect of proton-rich nuclei with N=Z is that neutrons and protons occupy 
the same shell-model orbitals. Consequently, due to the large spatial overlaps between 
neutron and proton single-particle wave functions, the proton-rich N=Z nuclei are expected 



to exhibit unique manifestations of proton-neutron (pn) pairing ||132| pL 36f| . 

At present, it is not clear what the specific experimental fingerprints of the pn pairing are, 
whether the pn correlations are strong enough to form a static pair condensate, and what are 
their main building blocks | 137|| . Most of our knowledge about nuclear pairing comes from 
nuclei with a sizable neutron excess where the isospin T—l neutron-neutron and proton- 
proton pairing dominate. Now, for the first time, there is an experimental opportunity to 
explore nuclear systems in the vicinity of the N=Z line which have many valence np pairs; 
that is, to probe the interplay between the like-particle and pn (T=0,1) pairing channels. 

This novel situation calls for the generalization of established theoretical models of nu- 
clear pairing. In spite of several early attempts to extend the independent quasi-particle 
formalism to incorporate the effect of pn correlations in light nuclei (see Ref. ||138|| for an 
early review), no symmetry- unrestricted calculations for np pairing, based on the isospin- 
projected quasi-particle theory, have been carried out. 

So far, the strongest evidence for enhanced pn correlations around the N=Z line comes 
from the measured binding energies. An additional binding (the so-called Wigner energy) 
found in these nuclei manifests itself as a spike in the isobaric mass parabola as a function 
of T 2 =|(iV — Z) (see the review [ |1 3Q| ] and Refs. quoted therein). The pn correlations are 
also expected to play a role in beta decay [|140| - |143l , deuteron transfer reactions ||144| , |145|| , 
structure of high spins ||138| , |146|J147|1 , and also in nuclear matter [|148j , |149| , j32] . 

The role of the T=0 part of the interaction on the presence of binding-energy irregularities 
near the N=Z has been recognized in Ref. 



152 . Recent calculations 151 



have revealed 

the rather complex mechanism responsible for the nuclear binding around the N=Z line. In 
particular, it has been found that the Wigner term cannot be solely explained in terms of 
correlations between the proton- neutron J—l, T=0 (deuteron-like) pairs (see Fig. [10]) . (For 



more discussion of this point, see also Ref. ||153|| .) 

Recently, the isospin structure of the density matrices and self-consistent mean fields has 



been discussed ||154j| in the HFB theory allowing for a consistent microscopic description of 
pairing correlations in all isospin channels. Theoretically, the pn pairing correlations have 
been studied by several authors in the HFB framework ||155| , |156|| . However, in Ref. ||154|| this 
has been done in the coordinate space allowing for the classification of generic proton-neutron 
mixing mean fields. The resulting HFB equations have interesting properties. For spherical 
nuclei, only the T=l and J=0 nucleonic pairs are allowed. (The presence of the T=0 and 
Jj^O pairs would necessarily lead to deformed mean fields.) Consequently, for the spherical 
symmetry, two cases can be considered. The first one corresponds to pn pairs coupled to 
S'=0 and is similar to the standard like-particle pairing. The second one, analogous to triplet 
pairing with S — 1, is more interesting because of breaking of the intrinsic parity. Here, 
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the HFB solution contains mixtures of proton and neutron states with different parities but 
with the same value of j, e.g., a mixing of the z/g 7/ / 2 and vr/7/2 orbitals. 
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FIG. 10. The strength of the Wigner term, W, extracted using binding energies calculated with 
the Ohoj shell model. Full shell-model calculations (filled circles) agree very well with experimental 
data (open circles). The results of shell-model calculations with the (J=l, T=0) two-body matrix 
elements removed ( J m ax=l variant, triangles) and with all T=0 matrix elements removed ( J m ax=7 
variant, diamonds), are also shown. The inset shows the values of W extracted from the ETFSI 
mass formula [|150|]. They are practically zero for all nuclei considered. (From Ref. [|151|.) 



B. Proton emitters 

Nuclei beyond the proton drip line are ground-state proton emitters. Initially, the parent 
nucleus is in a quasistationary state, and the proton emission may be considered as a pro- 
cess where the proton tunnels through the potential barrier. In most cases, the combined 
Coulomb and centrifugal potentials give rise to barriers which are as large as ~15 MeV. 
Consequently, the associated lifetimes, ranging from 10~ 6 sec to a few seconds, are suffi- 
ciently long to obtain a wealth of spectroscopic information. Experimentally, a number of 
proton emitters have now been discovered in the mass regions v4~110, 150, and 170 (see 



Refs. ||157| , |158|| and references quoted therein). It is anticipated that new regions of proton- 
unstable nuclei will be explored in the near future using radioactive nuclear beams. 

The width of the proton resonance can be estimated through the distorted-wave Born 
approximation (DWBA) | |159| ]. From the decay width one can obtain the half-life of the 
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proton emission, t\n. The proton resonances are extremely narrow, Y ~ 10 22 — 10 15 MeV, 



hence it is difficult to calculate their widths directly. (It is worth noting that in Ref. ||160 
an attempt was made to calculate the proton emission width by solving the Schrodinger 
equation in the complex plane.) 

Recently, half-lives of the spherical proton emitters have been calculated in Ref. ||161| 
in the core-plus-proton approach assuming Va p to be a sum of a simple nuclear optical 
Woods-Saxon potential and the Coulomb potential. Three different methods have been 
used: DWBA, the modified two-potential approach of Gurvitz ||162|| , and the semi-classical 
approximation WKB. (Because of its simplicity, the WKB approach has been widely used 



to study spherical proton emitters ||163|| .) 



After computing the barrier penetration factor, the experimental spectroscopic factors 
Sp Xp could be determined as ratios of calculated and measured half-lives [|158|| . Theoretically, 
the spectroscopic factor measures the fragmentation of a single-particle orbital (n£j). In the 
BCS theory, it is given by w 2 , i.e., the probability that the spherical orbital (n£j) is empty 
in the daughter nucleus. 



1 




FIG. 11. Correlation between the experimental proton spectroscopic factors S^ p deduced from 
measured ground-state proton emission half-lives of odd-Z, even-N proton emitters and theoretical 
values obtained in the BCS theory. The nuclei expected to be deformed are indicated by stars. 
(From Ref. pl| .) 

The correlation between experimental and theoretical spectroscopic factors obtained in 
the DWBA calculations is shown in Fig. |11] for odd-Z, even-iV ground-state proton emitters. 
The agreement between experiment and theory is good. For the two cases indicated by stars, 
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109 I and 113 Cs, the experimental values fall well below theoretical predictions. This suggests a 
strong fragmentation of the single-particle strength and/or increased tunneling probability as 
compared to spherical predictions. Indeed, both 109 I and 113 Cs are predicted to be deformed, 
and the "anomalous" proton half-life of 109 I has been reproduced by deformed calculations 



of Ref. fllEI 



In general, proton emission half-lives depend mainly on the proton separation energy 
and orbital angular momentum, but rather weakly on the details of intrinsic structure of 
proton emitters, e.g., on the parameters of the proton potential at least at a qualitative level 
(factors of 2-3). The weak sensitivity of ti/2 to the details of the optical proton potential 
has been discussed in Ref. |46| in the context of two-proton radioactivity. It has been shown 
that more than 94-99% of the WKB exponent comes from the region r>rg, which is almost 
solely determined by the combined Coulomb and centrifugal potentials. This suggests that 
the lifetimes of deformed proton emitters will provide direct information on the angular 
momentum content of the associated Nilsson state, and hence, indirectly on the nuclear 
shape. 

Proton radioactivity is an excellent example of the elementary three-dimensional 
quantum-mechanical tunneling. Experimental and theoretical investigations of proton emit- 
ters (or theoretically predicted ground-state di-proton emitters) will open up a wealth of 
exciting physics associated with the residual interaction coupling between bound states and 
extremely narrow resonances in the region of very low density of single-particle levels. 



V. CONCLUSIONS 

An experimental excursion into new territories of the chart of the nuclides will offer many 
excellent opportunities for traditional nuclear structure. This may include new regions of 
quadrupole and octupole deformation, new regions of shape isomers, including superde- 
formations, new combinations of magic or semi-magic closures, and many others. For a 



comprehensive review of these possibilities, the reader is referred to Ref. [165]. 

In trying to see the phenomena of a "new physics", we should ask the fundamental 
question of "how far is far"? Experiments with radioactive beams are going to be long 
and difficult, and many examples of nuclear exotica discussed in this paper (especially those 
concerning very neutron-rich systems) are clearly out of reach, even assuming most optimistic 
experimental scenarios. The hope is, however, that some of the effects associated with the 
loose binding will be seen as deviations from smooth systematic trends or will show 
up at higher excitation energies closer to the particle threshold, as in the example of the 
analog states | |12(j| , |i27| |. Theoretically, we are bound to adopt the strategy of going to the 
extreme values of N/Z in order to identify the qualitatively new phenomena, and then back 
down to experimentally achievable regions to see whether these phenomena can actually be 
observed. There is very little doubt that we are on the verge of the most fascinating fishing 
expedition; a lot of exciting physics will probably be caught already at the beginning of this 
journey. 

The main objective of this study was to discuss various theoretical facets of nuclear 
structure with radioactive beams. In particular, the unusual conditions created by the weak 
binding and the importance of the coupling to the particle continuum have been emphasized. 
The theoretical formalism has been applied to experimental observables; i.e., energy spectra, 
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masses, radii, surface thickness, and pair transfer form factors. It is demonstrated that these 
observables carry invaluable information that can pin down many basic questions regarding 
the effective nucleon-nucleon interaction. 

The analysis presented in this paper should be viewed as a useful starting point for future 
investigations. One of them is the coupling between vibrational and rotational modes and 
pairing fields in weakly bound nuclei. Another interesting avenue of exploration is the role 
of dynamics; e.g., the importance of the particle number conservation and isospin mixing. 
A fascinating and difficult research program is the microscopic description of excited states, 
especially those lying above the particle emission threshold. We are only beginning to explore 
many unusual aspects of the nuclear many-body problem offered by systems with extreme 
N/Z ratios. 
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